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ABSTRACT 


The  XYZ  Potential  Flow  Program  has  been 
modified  to  provide  an  improved  calculation 
of  streamlines  on  the  body  surface.  The 
improved  program  also  computes  the  geodesic 
curvatures  of  both  the  streamlines  and  the 
equipotential  lines  to  provide  input  for 
three-dimensional  boundary  layer  programs 
using  the  small  cross  flow  approximation. 


ADMINISTRATIVE  INFORMATION 


This  work  was  sponsored  by  the  Naval  Ship  Systems 
Command  under  Subproject  SF  14532106,  Task  15325,  Work  Unit 
1843-350. 


I.  INTRODUCTION 


The  value  of  boundary  layer  calculations  for  problems  of 
two-dimensional  or  axisymmetric  flow  has  been  amply  demonstrated. 
It  is  generally  believed  that  three-dimensional  boundary  layer 
calculations  will  be  of  similar  value.  A number  of  methods  have 
recently  been  developed  for  calculating  three-dimensional 
boundary  layers  with  small  cross  flow.  These  methods  require 
that  the  values  of  the  velocity  ard  of  the  geodesic  curvatures 
of  both  the  equipotential  lines  and  the  streamlines  be  specified 
along  streamlines  on  the  body  surface. 

The  XYZ  Potential  Flow  Program  has  recently  been  modified 
to  provide  input  data  for  three-dimensional  boundary  layer 
calculations.  These  modifications  are  part  of  a larger  effort 
to  improve  the  XYZ  Potential  Flow  Program  by  including  curved 
surface  elements  and  linear  source  terms.  A hybrid  program 
including  parts  of  the  older  and  newer  versions  is  being  made 
available  now  because  of  the  great  interest  shown  in  three- 
dimensional  boundary  layer  calculations.  This  report  »a  designed 
to  be  used  with  NSRDC  Report  3892*  on  the  XYZ  Potential  Flow 
Program  and  is  not  self-contained. 

Perhaps  it  would  be  wise  to  emphasize  here  that  small 
cross  flow  boundary  layer  calculations  are  still  of  an  experimental 
nature.  For  example,  the  starting  conditions  and  the  point  of 
transition  from  laminar  to  turbulent  flow  must  be  correctly 
specified  if  accurate  results  are  to  be  obtained.  These 
conditions  are  generally  not  known  and  few  people  have  much 
experience  at  guessing  them.  Thus, this  report  is  directed  more 
to  those  interested  in  developing  boundary  layer  theory  than  to 
those  doing  design  studies. 


1 Dawson,  Charles  W.  and  Dean,  Janet  S. , "The  XYZ  Potential  Flow 
Program,"  NSRDC  Report  3892,  June  1972. 
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II.  THREE-DIMENSIONAL  BOUNDARY  LAYER  METHODS 


Figure  1 shows  a three-dimensional  body  with  two  stream- 
lines. Note  that  the  streamlines  converge  or  diverge  and  curve 
as  they  pass  along  the  body  side.  The  convergence  or  divergence 
of  the  streamlines  can  be  modeled  for  each  streamline  by  an 
axisymmetric  body  as  shown  in  Figure  2.  In  other  words,  when 
the  radius  of  the  axisymmetric  body  changes  in  the  right  way, 
the  streamlines  will  converge  or  diverge  at  the  same  rate  as 
they  do  at  the  particular  streamline  on  the  three-dimensional 
body.  The  rate  of  convergence  will  change  from  streamline  to 
streamline,  so  a different  body  is  needed  for  each  streamline. 

Thus  as  a first  approximation  the  boundary  layer  on  the  three- 
dimensional  body  may  be  obtained  by  solving  a set  of  axisymmetric 
boundary  layer  problems. 

The  curvature  of  the  streamlines  on  the  body  surface 
indicates  the  presence  of  a pressure  gradient  across  the  stream- 
lines so  that  there  will  be  a cross  flow  within  the  boundary 
layer.  When  the  usual  small  cross  flow  approximation  is  used, 
there  is  an  equation  for  the  streamwise  boundary  layers  and  an 
equation  for  the  cross  flow.  The  cross  flow  depends  upon  the 
streamwise  flow  but  the  streamwise  flow  is  usually  treated  as 
being  independent  of  the  cross  flow.  Thus,  the  streamwise  flow 
computed  from  the  small  cross  flow  approximation  will  be  the 
same  as  that  computed  by  an  axisymmetric  program.  A computer 

program  of  this  type  developed  at  Douglas  Aircraft  Company  by 

. . . 2 

Cebeci,  Mosmskis  and  Kaups  is  available  at  NSRDC. 


^ Cebeci,  T;  Mosinskis,  G. J. ; Kaups,  K,  "A  General  Method  for 
Calculating  Three-Dimensional  Incompressible  Laminar  and 
Turbulent  Boundary  Layers  I.  Swept  Infinite  Cylinders  and 
Small  Cross  Flow,"  Douglas  Aircraft  Company  Report  No. 

MDC  J5694  , Nov  1972. 
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Figure  1 - Streamlines  on  a Three-dimensional  Body 


Figure  2 - Axisymmetric  Pipe  Model  of  a Three-dimensional 
Body  for  One  Streamline 


III.  NATURE  OF  THE  PROBLEM 


The  XYZ  Potential  Flow  Program  computes  an  approximate 
solution  to  the  problem  of  three-dimensional  potential  flow 
over  an  arbitrarily  shaped  body.  The  body  surface  is 
represented  by  a set  of  plane  quadrilaterals.  In  the  solution 
the  values  of  the  velocity  are  computed  at  the  centroids  of  the 
quadrilaterals.  A subprogram  which  computes  the  on-body  stream- 
lines is  given  only  the  quadrilaterals  and  the  velocity  at  the 
centroids.  The  earlier  subprogram  assumed  that  the  velocity 
was  constant  in  each  of  the  quadrilaterals.  The  resulting 
streamline  data  was  quite  rough  and  did  not  include  any 
information  on  the  curvature  of  the  streamlines  or  equipotential 
lines.  The  new  subprogram  uses  data  from  two  of  the  neighboring 
quadrilaterals  and  thus  can  include  linear  terms  in  the 
approximation  to  the  velocity.  As  a result  the  streamline  data 
is  much  smoother  and  the  geodesic  curvatures  of  the  streamlines 
and  equipotential  lines  can  be  computed.  Using  the  geodesic 
curvatures  the  metric,  or  effective  radius  of  the  body,  along 
the  streamline  is  computed.  This  data  can  then  be  used  with 
either  a small  cross  flow  boundary  layer  program  or,  ignoring 
cross  flow,  with  an  axisymmetric  boundary  layer  program. 
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IV.  RESTRICTIONS  ON  THE  INPUT 

The  input  is  the  same  as  that  described  in  NSRDC  Report  3892. 
However,  there  are  four  restrictions.  These  restrictions  are 
required  because  the  quadrilaterals  must  be  organized  so  that 
neighboring  quadrilaterals  are  readily  available.  The  quadri- 
laterals are  organized  in  groups  of  four  about  a common  corner 
point  as  in  Figure  3.  Thus  there  must  be  an  even  number  of 
quadrilaterals  in  both  the  N and  M directions  in  each  of  the 
sections  of  the  body. 

The  sides  between  quadrilaterals  are  used  to  define  the 
quadrilateral  coordinate  systems  and  serve  as  the  axis  of  rotation 
when  the  surface  is  flattened  to  facilitate  the  numerical 
differentiation  of  the  velocity.  Thus,  the  common  corner  point 
of  a group  of  four  quadrilaterals  must  not  coincide  with  any 
other  corner  point.  The  third  restriction  is  that  each  set  of 
four  quadrilaterals  must  have  at  least  seven  distinct  corner 
points  to  allow  a curved  surface  to  be  fitted  to  the  points. 


Figure  3 - Ordering  of  the  Quadrilaterals 
For  A Right  Handed  Coordinate  System 


Thus  only  two  of  the  four  quadrilaterals  may  degenerate  into 
triangles  by  having  two  of  their  corner  points  coincide. 

The  last  restriction  also  is  a result  of  the  way  the 
curved  surface  is  fitted  to  the  points.  The  angle  between 
the  normal  vectors  of  two  quadrilaterals  in  a group  of  four 
must  be  less  than  90°  and  perferably  less  than  45°.  If  a 
sharp  edge  is  required  the  input  should  be  arranged  so  that 
the  edge  is  along  the  boundary  of  the  groups  of  four  and  not 
through  the  center. 


V,  COMPUTATION  OF  THE  PARTIAL  DERIVATIVES  OF  THE  VELOCITY 


The  quadrilateral  is  placed  in  a standard  position  with 
an  upper  neighbor  and  a right  neighbor  as  shown  in  Figure  2. 
The  other  neighbors  are  ignored.  The  quadrilateral  coordinate 
system  is  designated  by  x^,  yq,  z^,  the  right  neighbor 
coordinate  system  is  designated  by  xr,  yr,  zr  and  the  upper 
neighbor  coordinate  system  by  xu,  yu,  zu.  Note  that  x^  is 
parallel  to  the  upper  boundary  of  the  quadrilateral  and  xr 
is  parallel  to  the  right  boundary  of  the  quadrilateral. 


Figure  4 - Standard  Position  for  a Quadrilateral 
and  its  Neighbors 

In  general,  the  neighbors  will  not  be  in  the  same  plane  as 
the  quadrilateral.  Thus,  the  first  task  is  to  compute 


appropriate  velocity  and  position  vectors  for  the  neighbors  in 
the  quadrilateral  coordinate  system  as  though  they  were  in  ' ne 
quadrilateral  plane.  This  is  done  as  follows: 


1.  The  cosine  of  the  angle  of  rotation  between  the 
quadrilateral  plane  and  the  plane  of  the  upper  neighbor  (cos  e^) 
is  found  from  the  dot  product  of  their  unit  normal  vectors. 

2.  The  upper  neighbor  velocity  vector  is  transformed  into 
the  quadrilateral  coordinate  system.  The  yq  component  is 
divided  by  cos  0U  to  correct  for  the  effect  of  surface 
curvature  and  the  normal  component  is  set  to  zero.  The  result 
is  as  though  the  velocity  vector  had  been  rotated  about  the 
boundary  between  the  quadrilateral  and  the  upper  neighbor  into 
the  plane  of  the  quadrilateral. 

3.  The  coordinates  of  the  centroid  of  the  upper  neighbor 
are  transformed  into  the  quadrilateral  coordinate  system.  An 
approximate  correction  for  the  effect  of  surface  curvature  is 
made  to  the  y coordinate. 

q 


1 2 2 i 

y corrected  = -r  [4  {y  + z ) + y + y /cos  0 ] 
-'q  6 q q ' q q u 


(1) 


This  equation  was  obtained  as  follows.  The  arc  length  along  a 
curve  in  the  y,  z plane  is 


i --  / 

^0 

When  approximated  by  Simpson’s  Rule,  f becomes 


Now  Yl-y0  = yq 

( J5-)  = o 

dy  y0 


= tan  9u  thus  |l+(^) 
ay  yi  u 


l/cos  0. 


u 
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is  not  known,  the  term 


Since  (|i)yj 


^rV  [ 


1+  (^)2y 


]* 


is  replaced  by  the  chord  length 


Thus  equation  (1)  is  obtained,  and  the  normal 
component  can  now  be  ignored . 


4.  The  cosine  of  the  rotation  angle  (cos  0 ) between 

r 

the  planes  of  the  quadrilateral  and  its  right  neighbor  is 
computed  from  the  dot  product  of  their  unit  normal  vectors. 


5.  The  right  neighbor  velocity  vector  is  transformed  into 
the  right  neighbor  coordinate  system  and  the  y^  component  is 
divided  by  cos  0^.  The  resulting  vector  is  transformed  into 
the  quadrilateral  coordinate  system  and  the  normal  component  is 
set  to  zero.  As  with  the  upper  neighbor  this  has  the  effect  of 
rotating  the  vector  about  the  right  boundary  into  the  quadri- 
lateral plane. 

b.  The  vector  from  the  quadrilateral  centroid  to  the  right 
neighbor  centroid  is  transformed  into  the  right  neighbor 
coordinate  system.  The  y^  component  is  corrected  for  surface 
curvature  and  for  the  transformation  to  the  quadrilateral 
coordinate  system. 

2 2 i 

y„  corrected  = [4(y  + z ) +y  +y  /cos  0 ]/S  cos  0 (2) 

Jr  rr  r r r r 

The  normal  component  is  set  to  zero  and  the  resulting  vector  is 
transformed  into  the  quadrilateral  coordinate  system. 


The  velocity  and  coordinates  are  now  known  for  three  points 
in  a plane.  Thus  the  two  tangential  components  of  the  velocity 
may  be  approximated  by 


U = Uq  + xq  + u2  yq 


V = Vq  + Vx  xq  ♦ V2  yq 


(3) 
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where  U is  the  xq  component  of  the  velocity 

V is  the  yq  component  of  the  velocity 

The  equations  from  the  three  points  readily  yield  the 
three  unknowns  in  each  approximation.  Thus  (U^,V^)  and 
(U2<V2>  which  are  the  approximate  partial  derivatives  of  the 
velocity  with  respect  to  xq  and  y^  are  obtained. 
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VI.  COMPUTATION  OF  POINTS  ON  A STREAMLINE 


The  streamline  computation  begins  at  a given  point  in  a 

specified  quadrilateral  with  the  velocity  vector  V = (U,V) 

given  by  equations  (3).  A stream  function  is  desired  which  will 

have  a constant  value  of  zero  along  the  streamline.  Since  the 

• • • 

divergence  of  V is  not  zero,  an  additional  function  p must 
be  found  such  that  the  divergence  of  p V is  zero.  Then  for 
a particle  following  a streamline: 

dy  _ dy/dt  _ V _ pV 
dx  dx/dt  U pU 

Thus  a new  vector  field  is  constructed  whose  streamlines  are 
identical  to  those  of  the  velocity  field  but  whose  divergence 
is  zero.  The  function  p is  given  by  the  series: 


u <Vv2) 

p(x  .y  ) = i - — * — -s — x 


q q 


-■* — y 
U^+V 

q q 


v (u ,+v,) 
a -1  ’2  y ♦ 

u2+v2  q 
q q 


(4) 


The  resulting  stream  function  is: 


SF<VV  = SF0  - VqXq+  Uqyq“ 


U^V(n(U.+V7) 

V - _a_s l l 


1 2 

x. 


2 T 
Uz  + v 
q q 


■f  (5) 


U V (U  +V.)" 

< 

U^CU^+V.) 

q q 1 2 

ia  + 

,,  q i 2 

i U +V 

7 

U1  2 2 

1 u +V 

L q q 

L q q J 

y 

q q 
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The  constant  value  SFq  is  chosen  so  that  the  stream  function 
is  zero  at  the  specified  point.  Since  the  velocity  is  known 
only  through  the  linear  terms  it  is  consistent  to  truncate 
the  stream  function  after  the  quadratic  terms. 

Values  of  SF  at  the  four  corner  points  of  the  quadri- 
lateral are  now  computed.  By  comparing  the  signs  of  SF  at 
adjacent  corners,  the  sides  through  which  the  streamline  passes 
are  determined.  If  the  streamline  passes  through  a side,  the 
value  of  SF  is  computed  at  the  midpoint  of  the  side.  From 
these  three  points  the  intersection  point  is  computed  from  a 
three  point  interpolation  formula. 

The  parameter  t(0<t<l)  is  defined  so  that: 


X = X1  ♦ t (Xg-X^) 

y = v1  + t (y3-y1> 


(6) 


Here  x,  y 

is 

the 

xl’yl 

is 

one 

X3’y3 

is 

the 

x2  ,y2 

is 

the 

intersection  point, 
corner  point. 

other  corner  point. 

middle  point. 


Then : 

0 = SF(x,y)  = 2 SF1  <t-i)(t-l)  - 4 SF2  t(t-l)  + 2SF3t(t-i)  (7) 

The  root  of  this  equation  is  chosen  so  that  0 < t < 1.  Note 
that  if  there  had  been  two  roots  between  0 and  1,  then  SF 
would  have  the  same  sign  at  both  corner  points  and  those 
intersection  points  would  be  ignored  completely. 

In  general,  two  intersections  will  be  found  for  the  entire 
quadrilateral,  one  where  the  streamline  enters  and  one  where  it 
leaves.  It  is  possible,  however,  for  there  to  be  four  inter- 
section points.  The  next  point  on  the  streamline  is  chosen  from 
the  intersection  points  as  follows.  A quantity  Q is  computed 
for  each  intersection  point  by  taking  the  dot  product  of  the 
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vector  from  the  starting  point  to  the  intersection  point  with  the 
velocity  vector  at  the  centroid  of  the  quadrilateral  and  then 
multiplying  it  by  the  direction  sign.  (The  direction  sign  is 
♦1  if  the  streamline  is  being  traced  downstream  and  -1  if  it  is 
being  traced  upstream.)  The  intersection  point  with  the  maximum 
positive  value  of  Q is  chosen.  However,  if  the  largest  Q 
is  less  than  or  equal  to  zero,  none  of  the  points  is  acceptable 
and  the  program  searches  for  another  quadrilateral  through 
which  the  streamline  might  pass. 

Since  each  of  the  points  representing  a streamline  is 
located  at  the  boundary  between  two  quadrilaterals,  the  values 
of  the  velocity  and  geodesic  curvatures  at  each  point  can  be 
computed  separately  for  the  two  quadrilaterals.  The  average 
of  the  two  values  at  a point  is  taken  as  the  value  at  that 
point . 

The  streamline  is  first  traced  in  the  downstream  direction 
and  then  in  the  upstream  direction.  After  the  entire  stream- 
line has  been  obtained  the  arc  length  along  the  streamline  is 
computed  so  that  it  starts  at  zero  at  the  upstream  end.  Points 
which  are  very  close  together  are  combined.  Such  points  occur 
when  the  streamline  just  cuts  across  a corner  of  a quadrilateral. 
If  the  distance  between  two  points  is  less  than  1/8  the 
distance  between  their  neighbor  in  front  and  their  neighbor  in 
back,  they  are  combined  and  average  values  of  their  velocity, 
position  and  curvatures  are  used  for  the  new  point. 
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VII.  COMPUTATION  OF  THE  GEODESIC  CURVATURES 
(Kj  AND  K2)  AND  THE  METRIC  COEFFICIENT  H2 


The  unit  vector  tangent  to  the  streamline  is 


iUF  jVt  kW 

(U2  + v2  + wV 


where  U , V , W 
1,  k 


are  the  components  of  the  velocity  vector. 

are  unit  vectors'  in  the  three  coordinate 
directions . 


The  curvature  of  the  streamline  is  defined  by 

K = d T/d l 


where  l is  the  arc  length  along  the  streamline.  The  geodesic 
curvature  is  the  component  of  K on  the  body  surface 

K = K • (T  x N) 

-V  g 

where  N is  the  unit  normal  vector  to  the  surface. 

For  a plane  quadrilateral  and  with  vectors  written  in  the 
quadrilateral  coordinate  system  we  have 


T = (i U + j V)  / (U2  + V2)i 

37  1 (u  h ' v 57)/<u2',v2>! 

K = (iV-jU)CU(VU  - U V )+V(VU  -UV  )]  / (U2+V2)2 

xx  y y 

where  subscripts  x and  y denote  partial  differentiation 

•*  7 7 i 

T x N = (-i  V + j U)  / (U  +V  ) 

Thus  K?»  the  geodesic  curvature  for  the  streamline  is  obtained. 

k2  = [ucuvx-  vux)+  V(UVy-  VUy)]/  (U2*V2 ) 2 (8 
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The  geodesic  curvature  for  the  equipotential  lines,  , is 
found  by  simply  replacing  U by  V and  V by  -U  in 
equation  8. 

3_ 

K,  = [-V(VU  - UV  )+  U(VU  - UV  )]  / (U2+V2)2  (9) 

i xx  y y 

The  metric  coefficient  is  computed  from  the  values 

of  K, . By  definition 
1 

-1  3H2 
ki  ' h1  h2  aT 

is  defined  to  be  one  along  the  streamline  so: 

9H2 

h2  ki  = " t r 

This  equation  is  approximated  by: 

H?U  + U ) - H 2U)  - [H  2U)  Kx(f)  + H2(f  ♦ U)  Kj_(t  + A*  ) ] 

or  2 - A£  K U ) 

H2(f+Af)  = H?(f)  j + Af  K 1(£  + Af ) (10) 

where  A t is  the  chord  length  of  the  segment  of  the  stream- 
line  passing  through  the  quadrilateral. 

Note  that  is  determined  only  to  within  a multiplication 

constant.  This  constant  was  arbitrarily  chosen  so  that  H2 
is  one  at  the  specified  starting  point  for  the  streamline. 

Thus,  for  an  axisymmetric  body  H2  will  be  proportional  to  the 
radius  but  will  not  i ->  general  equal  the  radius. 


VIII.  POSSIBLE  DIFFICULTIES 


There  are  several  possible  difficulties  in  the  computation 
of  a streamline  that  can  cause  the  program  to  end  the  stream- 
line in  the  wrong  place.  These  are: 

1 . When  the  projection  of  the  starting  point  onto  the 
piano  of  the  starting  quadrilateral  lies  outside  the  starting 
quadrilateral,  there  are  two  possiblities  which  may  occur. 

Tho  streamline  will  be  traced  in  only  one  direction  if  it 
passes  through  the  starting  quadrilateral  or  it  will  not  be 
traced  at  all  if  it  does  not  pass  through  the  starting  quadri- 
lateral . 

2.  If  the  streamlines  are  converging  from  both  sides 
toward  the  boundary  line  between  two  quadrilaterals , it 

may  happen  that  due  to  numerical  approximations  the  computed 
streamline  will  intersect  the  boundary.  Then  the  program 
cannot  find  another  quadrilateral  in  which  to  continue  the 
streamline  and  will  end  the  streamline  at  that  point. 

3.  The  present  program  does  not  provide  for  the 
continuation  of  a streamline  across  a plane  of  symmetry.  Hence, 
when  a streamline  intersects  a plane  of  symmetry  it  will  be 
terminated . 


IX.  OUTPUT 


The  edited  output  includes  the  following  information: 

(See  Appendix  B for  a sample  of  edited  output.) 

1.  The  problem  identification  and  information  from  the 
parameter  card. 

2.  Information  concerning  each  of  the  quadrilaterals 
specified  by  the  input.  This  includes 

a.  Warning  messages  about  possible  errors.  A 
message  "Questionable  Point  - Poor  Fit"  replaces 
the  "LARGE  D"  message  described  on  page  39  of 
NSRDC  Report  3892. 

b.  M,N  - the  indices  for  the  first  corner  point. 

c.  P - the  quadrilateral  number  in  the  total  array  of 
quadrilaterals . 

d.  XI,  Y1 , Z1  - the  coordinates  given  for  the  first 
corner  point.  (Point  M,N.) 

e.  X2,  Y2,  Z2  - the  coordinates  given  for  the  second 
corner  point.  (Point  M+1,N.) 

f.  X3,  Y3,  Z3  - the  coordinates  given  for  the  third 
corner  point.  (Point  M+1,N+1.) 

g.  X4  , Y4  , Z4  - the  coordinates  given  for  the  fourth 
corner  point.  (Point  M,N+1.) 

h.  XF,YP,ZP  - the  coordinates  computed  for  the  centroid. 

i.  XN  YN,ZN  - the  components  of  the  normal  vector. 

j.  A - the  area  of  the  quadrilateral. 

k.  FL  - xhe  maximum  distance  of  a corner  point  from 
the  centroid  of  the  quadrilateral. 
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1.  CZ1 ,CZ4 ,CZ5 , CZ6  - the  coefficients  of  a local 

quadratic  fit  of  the  body  surface.  These  are  not 
used  in  the  present  version  of  the  program  except 
for  checking  input. 

3.  Information  about  the  convergence  of  the  iterations  for 
computing  the  source  density.  This  information  includes 

a.  The  sum  of  the  absolute  values  of  the  changes  in 
the  source  density  from  the  last  iteration. 

b.  A,  B^ , and  - the  extrapolation  coefficients 
computed  from  the  last  five  iterations.  (Once 
every  five  iterations.) 

c.  A message  indicating  extrapolation  has  been 
performed . 

4.  The  edit  of  the  final  solution  includes,  for  each  of 
the  three  basic  flows 

a.  The  point  number  - the  same  number  as  the  P in 
Part  2.c  of  the  output. 

b.  X,Y,Z  - the  coordinates  of  the  centroid. 

c.  VX,VY,VZ  - the  components  of  the  velocity  at  the 
centroid . 

d.  ABS  • V - the  absolute  value  of  the  velocity. 

e.  CP  - the  pressure  coefficient,  CP  = 1 - V^/V^  . 

f.  The  source  density. 

g.  The  normal  component  of  the  velocity. 

5.  The  edit  of  the  solution  for  each  additional  flow 
includes  the  same  items  as  for  the  three  basic  flows  except 
that  items  f.  and  g.  are  omitted. 


19 


6.  The  output  for  each  on-body  streamline  includes  the 
following  quantities  for  each  of  a set  of  points  on  the  stream- 
line : 

a.  X,Y,Z  - the  coordinate  of  the  point. 

b.  VX , VY , VZ  - the  components  of  the  velocity. 

c.  CP  - the  pressure  coefficient. 

d.  K1,K2  - the  geodesic  curvatures. 

e.  H2  - the  metric  coefficient. 

f.  L - the  arc  length  along  the  streamline  from  the 
upstream  end  to  the  current  point  (X,Y,Z). 

g.  V - the  absolute  value  of  the  velocity. 


X.  USE  OF  THE  OUTPUT  FOR  BOUNDARY  LAYER  CALCULATIONS 


Most  boundary  layer  programs  require  input  consisting  only 
of  the  values  of  K2 , H2 , L,  and  V which  are  calculated  by  the 
XYZ  Potential  Flow  Program  and  described  in  Chapter  VIII.  It 
might  seem  then  that  using  this  potential  flow  data  to  do  a 
boundary  layer  calculation  would  be  very  simple.  Unfortunately 
there  are  several  points  which  must  be  considered  before 
proceeding  with  the  boundary  layer  calculation. 

First,  the  data  as  it  comes  from  the  XYZ  Potential  Flow 
Program  may  be  rather  rough  or  bumpy  and  may  not  contain 
enough  points  for  an  accurate  boundary  layer  calculation.  Thus, 
it  may  be  necessary  to  use  some  curve  fitting  routine  to  fit  a 
smoother  curve  through  the  data  and  thereby  increase  the  number 
of  data  points. 

Another  problem  is  chat  of  specifying  the  initial  velocity 
profile  and  the  point  of  transition  from  laminar  to  turbulent 
flow.  It  should  not  be  assumed  a priori  that  a method  which  is 
satisfactory  for  two-dimensional  or  axisymmetric  problems  will 
be  equally  satisfactory  for  small  cross  flow  problems. 
Experiments  with  several  methods  may  be  necessary  to  obtain 
satisfactory  results. 

Axisymmetric  boundary  layer  programs  require  the  radius  of 
the  body  as  input.  Thus,  when  an  axisymmetric  boundary  layer 
program  is  used,  the  radius  of  the  analogous  axisymmetric  body 
must  be  computed.  This  radius  will  be  a constant  times  the 
computed  by  the  XYZ  program.  There  are  two  points  to  consider 
in  choosing  this  constant.  First,  the  change  in  the  radius 
(AR)  between  any  two  points  on  the  streamline  must  be  less  than 
the  arc  length  (Af)  between  them.  Secondly, if  the  initial 
velocity  profile  is  to  be  determined  by  the  wedge  angle  at  the 
leading  end  of  the  streamline,  then  AR  for  the  first  two 
points  must  satisfy  the  equation 

AR  = a£  sin  0,  (11) 
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where  0 is  the  half  angle  of  the  wedge.  Thus  the  values  of 
R along  the  streamline  are  determined  for  the  entire  stream- 
line. Note  that  the  values  of  R which  are  determined  by 
equation  11  might  not  satisfy  the  first  condition  that 
AR  £ A l everywhere  on  the  streamline.  This  may  or  may  not 
cause  trouble, depending  upon  the  way  the  boundary  layer  program 
is  coded. 


APPENDIX  A 

DECK  SETUP  FOR  STANDARD  RUN  OF  SAMPLE  PROBLEM  SPHERE 


JOBNAME.CMiCOQOO,  T 1 0 3 , P ? • 
CHARGE.  CXXX,PP°PPPPPPP,Rr,,U. 
RFL, 55000. 

LIMIT, 353L. 

REQUEST  . T ApEL',,*DF. 

REQUEST, TAPFC4,*PF. 

RF  QUL  6 T , T A 0 5 7 , *®F  # 

!TriCH(»^l,C'<yvr>ri7<4)Arr,MI7  = l) 

irTACH(PF?,CXXXO"2V4TATr,Mi-*=l> 

ATTACH(OC7,CXXX3F3VnATir,MW=l) 

ATTACH(">c4,CXXXPF4\/nATc,MPil) 

ATTACH(PF7,CXXX3F7V40ATC,HR=1) 
ATTACH, MAT  «CLI  1AMMAT4. 

L 0 AO  ( M A T I 
PF1.  . 

RFL.USOOJ. 

PF?. 

RFL, 45033. 

PF3. 


* USr  RS  JOB  CAPO  • 

* USERS  CHARGF  CARO  * 

* CXXX  IS  USERS  10  • 

* pppppppppp  is  'jsrps  • 

* JOB  OROFR  MUM V P * 
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• COMPILATIONS  * 


• OMIT  IF  MATINS  WAS  * 

* COMPILFO  WITH  PF1  • 


RFL.F!  JJO. 

PF  4 . 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1 DTNSRDC  REPORTS.  A FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH 
NICAL  VALUE  THEY  CARRY  A CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES,  CONTAIN  INFORMATION  OF  A PRELIM- 
INARY. TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A CASE  BY  CASE 
BASIS 


